


o 

in 



> 

en 

o 
o 

0^ 






Relativistic two-body system 
in (l+l)-dimensions 



Norman Dombey 



Center for Theoretical Physics, University of Sussex, 
Palmer, Brighton BNl 9QH, UK 



and 



Fuad M. Saradzhev 



Institute of Physics, Academy of Sciences of Azerbaijan, 
H.Javidpr. 33, 370143 Baku, Azerbaijan 



Abstract 



The relativistic two-body system in (1 + 1) -dimensional quantum electrodynamics is studied . 

r^ It is proved that the eigenvalue problem for the two-body Hamiltonian without the self-interaction 

^ terms reduces to the problem of solving an one-dimensional stationary Schrodinger type equation 

^ ivith an energy-dependent effective potential which includes the 5-functional and inverted oscillator 

■ ■ fjarts. The conditions determining the metastable energy spectrum are derived, and the energies and 

widths of the metastable levels are estimated in the limit of large particle masses. The effects of the 

self-interaction are discussed. 



1 INTRODUCTION 

In the study of two-body systems in quantum theory we often use single-particle equations. For 
hydrogenlike systems, for instance, we assume that one of the particles is much heavier (proton) and 
then reduce the two-body problem to the problem of motion of the lighter particle (electron) in an 
external field of the heavier one. To get more exact solution we need to take into account two-body 
effects and start with a two-body equation. 

As shown in im, in (1 -f- l)-dimensions the single-particle Dirac equation allows no hydrogen. The 
equation has solutions for a continuous set of energies, and the probability of finding the electron 
infinitely separated from the proton remains finite at all times. Despite the attractive force, the 
electron and proton are not confined in a hydrogenlike system with discrete energy levels. That 
happens not only for hydrogen atoms with an infinitely heavy source of potential, but also for 
positroniumlike systems. 

In the present paper, we want to clarify to what extent the two-body effects influence the result of 
imi. We aim to study a relativistic two-body system in (1 + l)-dimensional quantum electrodynamics 
(QED) by making use of a two-body Dirac equation. Models in (1 -|- l)-dimensions are known to 
be useful as simpler models for discussion of many-body aspects of particle physics, in particular, 
spontaneous positron production by supercritical potentials 0, |^. 

To describe two-body systems we usually introduce a composite field, and there are two ways of 
deriving equations on this field. If we rewrite the action of the two-body system entirely in terms of 
the composite field, then we can require the action to be stationary with respect to the variations 
of this field only. This way leads to a single two-body equation @] - 0. However, if we first vary 
the action with respect to the individual fields, then we come to a pair of coupled equations on the 
composite field. The pair of Dirac equations formulation of the two-body problem was given in [0, ^ 
in the framework of the constrait approach. The main difference between the two ways is in the role 
of the relative energy ( or its conjugate variable, the relative time). While in the first way the relative 
energy drops out of the two-body equation automatically, in the pair of Dirac equations formulation 
it is eliminated by using the compatibility condition of the two equations and a special choice of the 
interaction potential. 

In our paper we follow the single two-body equation formulation and work in the first-quantized 
version of QED when both matter and electromagnetic fields are not quantized. We consider a system 
of two massive Dirac fields minimally coupled to a ^(1) or electromagnetic field. The electromagnetic 
field has no separate local degrees of freedom and can be eliminated between the coupled Maxwell- 
Dirac equations, but then nonlinear self-field terms must be included. In Sect. 2, we derive a 
relativictic two-body equation in the self-field QED^^^ defined on the line and give the Hamiltonian 
form of this equation. In Sect. 3, we find the eigenf unctions and the spectrum of the two-body 
Hamiltonian. We study in detail two cases: i) free motion; ii) the Coulomb interaction, and discuss 
the effects of the self-interaction. Sect. 4 contains our conclusions. 

2 TWO-BODY EQUATION 

For our system, the action is 

roc roo 

W= / dt dxC{x,t), (2.1) 



2 _ _ 1 

k=i ^ 

where {jj.^u = 0, 1), 7° = —icr2, 7^7"'^ = 7^ = 0-3, ctj {i = 1,3) are Pauli matrices. The fields ipk are 
two- component Dirac spinors, and tpk = '4'kl^- The partial derivatives are defined as Oq = d/cdt, 
di = d/dx. 

The electromagnetic field equations deduced from the action ( |2.1|) are 



d^F''^' = J'^, (2.2) 

where the total matter current 

2 

k=l 

is conserved, 9^J^ = 0. 

In the Coulomb gauge Ai{x,t) = 0, the equations (^.21) take the form 



dlAo = -J°, 
didoAo = J\ 

These two equations reduce in fact to each other because of the total current conservation and are 
solved by 

POO 

Aoix,t) = - dyDix,y)J%y,t), (2.3) 



where the Green's function is 

F>{x,y) = -\x-y\. 

In contrast with the situation on the circle, the electromagnetic field on the line has not a global 
physical degree of freedom and can be therefore eliminated from the action completely. If we insert 
into Eq.( ^.lD , we obtain the action in the Coulomb gauge as 

2 

dt dxY, ^kil^i^cd^ - mkc'^)iPk+ 

-00 J— CO 1 1 



fc=l 



1 /"OO roo roo 

+ -/ dt dx dyJ\x,t)D{x,y)J\y,t). (2.4) 

2 J —00 J —00 J ~oo 



We could vary the action ( |2.4| ) with respect to individual fields ipi and ip2 separately. This 
results in non-linear coupled Hartree-type equations for these fields. Instead, we use a relativistic 
configuration space formalism |^, Q to take into account the long-range quantum correlations. We 
define a composite field $ by 

$(Xi,t|x2,t) = ^l(Xi,t) <S)ij2{x2,t) 

which is a four-component spinor field. The configuration space (xi, 0:2) is two-dimensional Euclidean 
space R^. 



We can rewrite our action ( |2.4| ) entirely in terms of the composite field $. The resultant action 

isfl 



W[<l>, A]= dt dxi / rfx2$(xi, t|x2, t){(c7>(i),^ - m^c^) ® 7° + 7° ® (c7>(2),m - ^20^) 



+ ^(7° ® 7°)(ei0S' + e20g^) + 6162(7° ® l'')D{x,, X2)}<l>ix^,t\x2, t), (2.5) 



where 



and 



P(^),^l = i^ 



d 



dx^ 



A" 



0^f/(x,t) = ei / rfy/ rfzD(x,2;)$(2;,t||/,t)(7°®7°)$(^,t|y,t), 

J —OQ J —OQ 

/OO /"OO _ 

dy / dzD{x, y)$(z, tly, t)(7° ® 7°)*(^, ib, t), 
-00 J —00 

the self-potentials 0(|\^ being non-linear integral expressions. The spin matrices are written here in 
the form of tensor products ® , the first factor always referring to the spin space of particle 1, the 
second to particle 2. 

Now we require the action ( |2.5| ) to be stationary not with respect to the variation of the individual 
fields but with respect to the composite field only. This leads to the following two-body wave equation 

{iY^^{l),^i - mic) (g) 7° + 7° (g) (7''7r(2),^ - m2c) 

+ ^i^(7° ® l')Dixu X2)Mx,,t\x2, t) = 0, (2.6) 

c 

where the generalized (kinetic) momenta vr(j)^^ are given by 

^ ^ I ^* /I self 

with 

/I self J,self /I self J>self 

^(1),0 = r(l) ; ^(2),0 = <P(2) ) 



and 



/I self /I self r\ 



In the center of mass and relative coordinates 



n = 7r(i) + 71(2) , TT = 7r(i) - 7r(2) , 

P = P{1) + P{2) , P = P{1) - P(2) , 

X_|_ ^ X\ ~T X2 J X^ ^ Xi — ^2; 



the function D(xi,X2) becomes 



D{xi,X2) = D_{x-) = -\x- 



i.e. depends only on the relative coordinate X- and is symmetric. 

The two-body equation , first without the self-field terms, takes the form 



r^P^ + Pp^ + ^1^(7° ® 7°)^- (a;-) - micJ ® 7° - m2C7'^ 



$(x_,t|x+,t) = 0, 



where 

P = ^(7^® 7° -70^7'^), 

and I is identity matrix. We see that kP vanishes which means that the relative energy po drops out 
of the equation automatically and we get 



popo _ pipi _ ^y ^ ^i^(^o ^ ^0)^_(^_) _ ^^^j ^ ^0 _ ^^^^0 



$(a;_,t|x+,t) = 0. (2.7) 



Thus we have only one time variable conjugate to the center of mass energy ip°, one degree of freedom 
for the center of mass momentum P^ and one degree of freedom for the relative momentum p^. Since 
P° is the " Hamiltonian" of the system, by multiplying ( p.7|) by Fq^^ we obtain the Hamiltonian form 
of the two-body equation 

Po$ = U+P^ + a-V^ - ^^D^ + f3imic + (32m2c) $, (2.8) 



with 

a± = -(alias), ai = 7^ ® /, as = / ® 7^ Pi=-f°®I, p2 = I^l^, 

and the relative and center of mass terms in the Hamiltonian Pq being additive: 

-P = Hem. + -"reh 

Hem. = a+P^, 

ifi-gi = a_p^ D_ + Piniic + [32m2C. 

c 

Eq. ( p.8|) has the form of a generalized Dirac equation, now a four-component wave equation. 
With the self-potential terms the Hamiltonian form of the two-body equation becomes 

Po$ = U+P^ + a.p" - -0_ - ^<pf^ - ^0^f + Amic + /Jamac) $, 

where 

0_ = eie2-D-. 

The self-potentials break in general the above mentioned additivity of the center of mass and relative 
parts of Pq. 

3 SPECTRUM 

Let us find the eigenfunctions and the spectrum of Pq. The equation for the eigenfunctions is 

{a+P^ + a_p^ + /?imic + /^zmac)^ = -{E + 0)$, (3.1) 

c 

where 

(9 (9 

pi = 2ih^^, p^ = 2ih- 



dx^ dx^ 



and 



+ ei0^f/ + 62^$;. 



If we denote the components of the composite field $ as 



then (|3.1|) reduces to the system of four equations: 

d 1 ~ 
2ih- — T]i (0 + E)r]i = -micr]3 - m2cr]2, 

OX-y. c 

d 1 ~ 

2ih- — r/4 + -(0 + E)t]4 = micr]2 + ^20773, (3.2a) 

OX+ c 

d I ~ 

2ih- — r]2 (0 + E)r]2 = -micr]^ - m2cr]i, 

OT- c 

d 1 ~ 

2ih- — r]3 + -(0 + E)r]3 = micr]i + 77120774. (3.2b) 

OX- c 

We see from these equations that 

7/^^(^,61,62) =774(^,61,62), (3.3a) 

r]2{E, 61, 62) = r]3{E, 61, 62), (3.3b) 

i.e. only half of all solutions of Eqs.(3.2a) - (3.2b) are independent and correspond to physical 
particles. 

It is more convenient to introduce the combinations 

V± = V2± Vs, 
X± = Vi±Vi, 

and use them instead of the original componenets rji [i = 1,4). The system of equations (3.2a-b) 

takes the form 

d 1 

2ih- — x+ - -fX- = {Am)cr]_, 

OXj^ c 

2ihj—x- - -fx+ = ~Mcr]+, (3.4a) 

ax+ c 

d 1 
2ih- — 77+ /7/_ = {Am)cx-, 

2^/1^— 77_ - ^fr]+ = -Mcx+, (3.4b) 

where M = mi + m2, Am = mi — m,2, and f = (p + E. 



while x+ ^"^^ V+ £''^6 related by 



where 



Without loss of generality, we can take masses equal to each other, mi = 7712 = m. Then the 
eigenfunctions X- ^"^^ V- ^-^^ determined by x+ ^^^ V+i respectively, 

T_?7+ = —2mcx+, (3.5a) 

T+X+ = -2mc77+, (3.5b) 

T± - -4^'^7^7^ - -/• (3-6) 

OX± J ox± c 

Acting on (3.5a) by T_ and on (3.5b) by T+, we easily decouple x+ and 1]+ and rewrite (3.5a-b) 
equivalently as 

T^T^T]^ = 4m^c^?7+, (3.7a) 

T^T+x+ = ^m^c\+. (3.7b) 

In what follows we assume for the eigenfunctions the ansatz 

r]±{x-,x+) = exp(-PcmX+ji]±{x^), 

X±{x_,x+) = exp(-PcmX+]x±{x-), 

i.e. separate their x+-dependent part from the x_-dependent one. Here Pcm is a momentum conjugate 
to the center of mass coordinate x+, so we can call it "the center of mass motion momentum". 

3.1 Free motion 

If we neglect both mutual and self-interactions, i.e. put = 0, then we get a system of two "free" 
particles. The relations (3.3a-b) become 

r]i{-E) = r]4{E), 
r]2{-E) = -r]s{E), 



i.e. the negative energy solutions of rji and 772 coincide correspondingly with the positive energy 
solutions of ri4 and 773. Therefore we may consider either positive and negative energy solutions of 771 
and r]2 or only positive energy solutions of all four equations (3.2a-b) as physical particles. 

For = 0, the operators T_ and T+ commute, so Eqs.(3.7a) and (3.7b) coincide with each other. 
Their solution is 

rj+ix) = sin ( —k,x_ 
\n 

, , Imc^E 



where 



E 



^2 _ 4c2p2 _ 4^2c4 



2c\\ E^-Ac^P^„ 
For the other two components, we have 



hi \h 



2CK /I 



cm 
X-^--) = -^2_4,2p2 ^-U"^- 



This solution exists for all energies for which 



/ P2 

\E\>2mc'^\ 1+ "' 



9 9 ' 



and there is a "forbidden" band in between (see Fig. la). 

3.2 Coulomb interaction 

Let us now consider our two-body system in the presence of the Coulomb interaction and neglect 
only self-field terms. The relations (3.3a-b) take the form 

T]i{-E,-eie2) = 774(^,6162), 
?72(--E, -6162) = -773(^,6162), 

i.e. the negative energy solutions of rji and ri2 coincide correspondingly with the positive energy 
solutions of ?74 and rj^ of opposite sign of 6162. 

For = 0_, the equation (3.7a) for 77+ reduces to the second order differential equation 

9^?7+ 1 f df \ dr]+ 1 2 "^^c^ p 



+ 772-;^/ ^+ - ^2-72 ^„2n2 ^+- (3-^ 



dxl f\dx^)dx_ ATi^c^' '+ Ti^ p-Ac^P^^ 
If we make in ( |3.8|) the substitution 

Ti+{x~) = v/-^(a;-), 

then we find for a the following Schrodinger type equation 



dx"^ 
with the "potential" 



+ V{x_)(j = Ka (3.9) 



,,, , 1 rfV 3 /I d/ V 1 ^2 4m2c^p2 



2/da;2 A\fdx_J Afi^c^'' h\P-4c'^P^^^ 
and the "energy" 



m'^c^ 



K= „ 
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The last term in the potential represents the center of mass motion contribution which vanishes for 
Pcm = as well as for all values of Pcm in the massless case. 
The explicit form of the potential for Pcm = is 



V{x^ 



-5(s_) +V{x. 



where s = ^^, 

6162' 



and 



V{x^] 



4(|x_| + s)2 16rc2 



(eie2)^(|x_| + s)^ 



(3.10) 



The potential V{x-) has several pecularities. First, it contains a 5- functional part with coefficient 
(— 1/s) , positive ( for 6162 > and E < or 6162 < and E > ) or negative ( for 6162 > and 
ii^ > or 6162 < and E < ). The form of V{x-) for different signs of the (5-function coefficient is 
shown in Figures 2 and 3. Secondly, its regular part V{x-) includes the inverted x^-potential (the 
last term in Eq. (|3.10| ) ). Such kind of potentials is known to appear in barrier penetration problems, 
splitting in double wells, and tunneling out of traps [|l^]- [Q. The systematic study of the inverted 
oscillator is given in |]13|. The presence of the inverted x^-potential makes V{x-) nonvanishing at 
|x_| — >■ cx) and indicates that the particles are not confined in a stable system. 

The regular part V{x^) is symmetric with respect to x_, V{—x^) = V{x. 
non-zero a:_, |a;_| <^ |s| , V(x_) is approximately linear 



For very small 



V(x_ 



Vn 



3 1 1 



(6162) 



x_ 



where 



Vo = V{x^ 



, 3 1 

^ 4s2 



6162 2 2 



The value Vq is positive for E'^ < ^hc\eie2\ and negative for E'^ > ^hc\eie2\- 

While cr(a;_) is continuous for all x_, its ffist derivative da/dx_ changes discontinuously at the 
point a;_ = 0. This is because of the 5- functional potential in the Schrodinger equation for a. If 
we integrate both parts of ( |3.9| ) over infinitely small interval (—£,£:), £ <^ 1, and then take the limit 
£ ^ 0+, we get the matching condition 



da , „, da 
+0)- 



dx 



dx. 



-0) = --a{0). 
s 



(3.11) 



The Schrodinger equation ( |3.9| ) taken without the center of mass motion contribution to the 
potential can be solved exactly. With the substitution 



aix. 



^'^'^-p(-'i^)^i')^ 



where z = |eie2(|a;_| + s)^, and away from the origin (a;_ j^ oi z j^ zq = \eie2s'^ = E^/eie2) the 
equation becomes 



d^u , z du i K 

dz"^ he dz he 6162 



'U = 0. 



The first independent solution of this equation is 



u^ = F{l + ip,2 



iz , 
he 



[3^ 



Khc 

6162' 



i.e. the confluent hypergeometric function. Tlie integral representation for ui is 0] 



Ml = 2 exp ( i—— ] Re 
2hcJ 



T{l + t(3) 



exp 



i'2Ei)0"'"^^'-^^--^^-Tj 



The definition and some useful formulae for the function G are given in Appendix. The asymptotic 
behaviour of the first solution is 



5 = aigT{l-il3). 



The second independent solution is 



«2 = -2exp(.^)lm 



exp 



('2Hj)©""^''('-"'--'^^i) 



Its asymptotic behaviour is 

0-2(1^1 -> OO) 

If we write the total solution 



T{l + i(3) 

|r(i + 2/?)| V2;ic '^ ^c y 

(J = Aai + B(72, 



where A and B are arbitrary constants, and use the matching condition ( ^.11| ) which in terms of z 

becomes 

a'izo) _ 1 1 

cr{zo) Azo' 

then we get the following relation between A and B : 

A Azoa'^izo) + a2izo) 



B 



4:Zoa[{zQ) + ai{zo)' 



the prime indicating the derivation with respect to z. 
Asymptotically the total solution behaves as 



a{\z\ 



OO 



hce 2^ 

W+m' 



-1/4, 



{{B - iA) exp [i—- + ipln— + iS 



2hc 
z6)}. 



he 



+ {B + iA) exp -i-— - i(3\n- 

V 2hc he / 

Regardless of values of the constants A and B, a does not vanish at \z\ ^ oo and can not represent 

bound states. For B + iA = 0, the solution behaves at infinity as a diverging wave. Such behaviour 

is specific for metastable states. The condition determining the metastable (or quasi-discrete) energy 

levels is then 

a[{zo) -ia'^jzo) ^ 1_ 

(Ji{zo) - ia2{zo) 4^0 ' 

With the expressions for ai and 0-2, we rewrite ( |3.12|) as 



(3.12) 



G(l-^/5,-z/?;^) 



2he zq 
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and then, using the relation between G and its derivative (see Appendix), finally come to 

,,„ /?(! - m G{2-i(3A-i(3-f^) _ 1 (3 

Metastable states are described by complex values of energy 

where Eq is the metastable level energy, while F is its width. For zq we get 

1 /rf r^ EqY 
6162 4 6162 

We can solve Eq.( p.l3| ) approximately for large values of m, m? ^ h\eie2\/c^- In this approxima- 
tion, \(3\ S> 1, so 

G{2 -1(3,1- 1(3- ^^ 



G(l-z/3,-z/3;^) 



1. 



Eq.( pJ.13D takes the form 



w-/5)4(|)^ + /5(i) 



which is solved by 



e1^^ ^ ±mcVi + Va- 



ra) 



1 

eie2\h 1 



J_n _ 5^3 . (6162^)' 

12^ 12 '^ mV 



^/3mc y/i + ys' 
and 

° '^ V6mc (1 + ^3)3/2' 
F(2)^v^mc2.(l + v^)3/2. 

There are therefore four metastable energy levels in the band between the positive and negative 
energy continuums, the first two at the energies Eq ^ ±mc^yl + y/S = ±1.65 x mc^ and the other 
two at the energies close to zero, Eq ^ 0± (see Fig. lb). For infinitely large values of mass, m -^ 00, 
the first two metastable energy levels turn into stable ones, F*^^) -^ 0, while the second two disappear, 
F(2) -. 00. 

For 6162 > 0, the positive energy metastable levels correspond to the relative motion in the 
potential V{x-) with s > 0, while the negative energy ones in the potential with s < 0. For 6162 < 0, 
the metastable energy levels are positive for s < and negative for s > 0. 

In the massless case (/3 = 0), the solutions cri and (72 become trigonometric functions 

0"! = 2hcz~^''^ sin ( — — ) , 
\2ncJ 

o"2 = 2hcz~^'^ cos ( — — 
\2hc 
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The analogue of the condition (|3.12|) is 



tg 



\2hcJ 



This condition does not fix Eq, while for the level width we get T = oo. This means that in the 
massless case there are neither discrete nor quasi-discrete energy levels and the energy spectrum is 
continuous. 

The spectrum for x+ can be derived analogously. It can be shown that the equation (3.7b) for 
X+ also reduces to the Schrodinger type equation with the potential 



U{x-) = V{x-) + 



Id'f fldf 



f dx^ \ f dx^ 
For Pcm = 0, the exphcit form of f/(x_) is 



4c^Pl + P 



2 o2 



c'P, 



df 



Ac^Pr. 



P P Ac^P^^-P\dx^ 



t/(x_) = — 5(2;_)+f/(x_), 
s 



ij{x. 



7 



-L /CiC2n2/| I I \2 



4(|x_| + s)2 16' he 
Acting along similar lines as above, we get the following metastable spectrum equation 






2nc 2 Zq 



For large values of m, this equation is solved by 



E, 



(1) 






1 . 1 ^ {eie2hy 

1 H — fi r ^ 

48^ 2V3' m^c^ 



and 



E, 



(2) 



2y/3mc 
6162^ 



i + Vs, 

1 

I + V3' 
1 

1 + v^ 



(3.14) 



2^/6mc 
r(2) ^ 2V2mc2 ■ 



The structure of the spectrum for x+ is the same as in the case of 77+. For infinitely large values of m, 
the spectrums for 77+ and x+ coincide exactly, while for large and finite values of m, the corrections 
of the order (1//5) and (l//9^) are different. 

3.3 Self-interaction 

The self-interaction makes the spectrum problem essentially more complicated. Let us give here 
a few comments concerning the effects of the self-field terms. 



12 



With the self-potentials 0/f/ , 0(2*/' ^^^ function / depends on both coordinates x_ and x^, and 
the operators T± acquire additional terms including the partial derivative (df/dx^). Moreover, we 
can not assume, as before in the study of the Coulomb interaction, that the center of mass motion 
is free, with a momentum Pcm- This results in infinitely large values of the self-potentials. 

Indeed, in terms of the components r]i (i = 1,4) the self-potentials take the form 

Co f°° f°° /I \ - ^ 

<Plf)ix) = y j_^ dy^ j^^ dy+D [x, -{y+ - yj)j J2'ntiy~^y+)vi{y-,y+)- 

For the free center of mass motion, 

Vi{y-^y+)vi{y~,y+) = Vi{y^)vi{y-), 

so the integrals /^ dy+D (x, |(i/+ ±1/-)) diverge. 

We can not assume also that the dependence of the components on the coordinates x_ and x^ 
is factorized, because for a general form of /(x_,x+) such factorization is simply not valid. Even in 
the case of the purely Coulomb interaction the factorization takes place only when the motion of the 
center of mass is free. To prove that let us use for a moment the following ansatz for the components 
r]+ and x+- 

V+ix-,x+) = ri+{x+)r]+{x_), 
X+{x~,x+) = x+{x+)x+{x-). 

With f = (p^ + E, the equation (3.5a) gives 

V+{x+) = X+{x+), 
while the equation (3.5b) becomes 

1 d'^x+ ^ Px+{x~) - 2/mcV(3;_) 
X+{x+) dxl -4a^c2x+(a;_) 

Since the left-hand side of this equation depends only on x+ and the right-hand one only on x_ , 
both sides must be equal to an arbitrary constant. Choosing the constant as {—Pcm/^^)i ^^ S^^ 

X+{.x+) = cxp (-Pc„,x+j , 

i.e. the factor corresponding to the free motion of the center of mass (if the constant is taken 
positive, say 1/-R^„, where Rem is a parameter of the dimension of length, then we come to the factors 
exp {±x+/Rcm) which diverge at positive or negative infinity and are therefore unacceptable ). 

The self-potentials are usually calculated by iteration procedure. To lowest order of iteration we 
solve the spectrum problem without the self-field terms. Then we substitute the solution obtained 
into the expressions for the self-potentials, calculate these potentials explicitly and use them in the 
next order of iteration. Thus, to get finite expressions for the self-potentials and to continue the 
iteration procedure we need at the lowest order, i.e. in the Coulomb interaction case, a general 
solution of the problem without the assumption of the free motion of the center of mass. 
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4 DISCUSSION 

1. We have shown that the spectrum problem for the two-body Hamiltonian in (l + l)-dimensional 
QED reduces to the problem of solving a system of two second-order partial differential equations . If 
the center of mass motion of the two-body system is assumed to be free, then these equations govern 
only the relative motion and take the form of one-dimensional stationary Schrodinger type equations 
with energy-dependent potentials which include the (5-functional and inverted ocsillator parts . 

We have solved the problem in the case of equal masses and the self-potentials neglected, and 
derived the conditions determining the metastable energy levels. We have estimated the energies and 
widths of the metastable levels for large values of mass. For the vanishing mass, neither stable nor 
metastable levels exist, and the energy spectrum is continuous. 

2. Our consideration on the basis of the two-body equation without the self-field terms does 
not change the result of the single-particle Dirac equation approach, namely, the nonexistence of 
hydrogenlike systems in (1 -|- l)-dimensions. However, the two-body equation, even with the self- 
interaction neglected, provides essentially new details: for limited times the particles can be confined 
in a metastable system characterized by quasi-discrete energy levels. For large values of the particle 
masses, the metastable system does not decay for a long time, and its spectrum is close to a discrete 
one. 

To treat the problem completely it is necessary to take into account in the two-body equation 
the self-potentials. It is also of a principal importance to consider the center of mass motion as a 
finite one, since for the free motion case the self-potentials take infinitely large values. This work is 
in progress. 
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APPENDIX 



The function G(q;,7; 2;) is defined as 

G{a, r, z) = ^AP^ / (1 + -y-t^-'e'dt, 
2m Jci z 

where the contour Ci comes from infinity (Ret -^ —00), goes round the point t = and then returns 
to infinity (see Fig. 4.). 

The asymptotic expansion of G for |2;| ^ 00 is 

z 2\z^ 

The relation between G and its derivative with respect to z is 

G'{a, 7; z) = 7rG{a + 1, 7 + 1; z). 

Z'^ 



Other relations are 



and 



G{a,^\z) = G{a + l,'y;z) G{a + 1,'y + 1; z), 

G{a,-f + l;z) = -G{a + l,-f + l;z) + G{a,r,z), 



G'{a, 7; z) = - [G{a, 7; z) - G{a + 1, 7; z)] , 
G'{a, 7; z) = - [G{a, 7 + 1; 2;) - G{a, 7; z)] . 
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Figure Captions 

Figure 1 

The spectrum of the two-body system for Pcm = 0; (a) free motion; (b) with the Coulomb 
interaction. The width of the metastable states is not shown. 

Figure 2 

The form of the potential V{x^) in the case of equal masses and without the self-field terms for 
s > 0, i.e. for 6162 > 0, -E > or 6162 < 0, £" < 0. Only the case E'^ < ^hc\eie2\ is shown. 

Figure 3 

The form of the potential V{x-) in the case of equal masses and without the self-field terms for 

2 



s < 0, i.e. for 6162 > 0, -E < or 6162 < 0, -E > 0. Only the case E'^ > ^hc\eie2\ is shown. 



Figure 4 

The contour Ci. 
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